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ABSTRACT 

Rotor unbalance is most common fault found in the rotating machines. Methods 
are adopted to analyze the position of unbalance and to bring its effect into acceptable 
limit. Vibration analysis is the most common technique used to analyze the rotor 
system. Research have been performed on rotor supported at both ends, however less 
study has been done for overhung rotor. In this paper the response of overhung rotor 
on isotropic support and anisotropic support subject to unbalance has been presented. 
Mathematical formulation is done using Finite Element method and equations are 
solved using MATLAB programming. The effect of unbalance has been studied on the 
bode plot. Forward and Reverse whirl are observed through Campbell diagram and 
mode shapes are plotted. 
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1. INTRODUCTION 


The unbalance response of the rotor system is an important aspect to consider while designing 
a rotor. Finite Element Method (FEM) has proved out to be of great importance in study of 
the dynamics of rotor system especially for unbalance response due to its accuracy in the 
results [1, 2], As the complexity of rotor problem increases the matrix size involved in the 
FEM analysis also increases and researches have been carried out for more than a decade to 
reduce the matrix size of the elements involves. Various Matrix reduction scheme have been 
proposed by researchers. Substructure synthesis method is one of the scheme which has 
gained importance which involves the division of the model into substructures followed by 
analyzing these substructures and finally compiling the results to give response of the overall 
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rotor system. An improved substructure synthesis method was presented for studying the 
unbalance response of multi-span rotor by introducing an exact matrix condensation scheme 
which reduces the size of the element matrix, consideration was made that a complicated rotor 
system consists mainly of multi-span rotors with less effect of coupling or bearings [3]. 

Model based identification techniques have been of keen interest of many researchers to 
study the dynamic behavior of the rotors. The most common technique for unbalance study is 
influence coefficient approach [4], The response of rotors has been studied combining two or 
more defects. Hertzian contact force, surface waviness, internal radial clearance etc. are some 
of the non linearities which have been combined with unbalance effect to study the response 
of rotor. Response maps, Frequency maps, orbit plots and Poincare maps are used to analyze 
the output response of the rotor system [5-8]. Dynamic behavior of the rotor gets changed by 
the presence of faults like unbalance. The changes are due to change in the equivalent loads 
that arise from moments and forces due to the damage. The effects of these equivalent loads 
are studied by applying them on the undamaged model. The model with these loads 
mathematically is known as fault model. Fault models of rotor can be used to compare the 
effect of other defect like unbalance which can help to determine the position and phase of the 
unbalance [9]. The method of equivalent loads has a limitation that the number of measured 
vibrations has to be kept less to avoid errors. This is achieved through theoretical fault method 
and vibration minimization method. These methods can measure unbalance fault even for two 
degree of freedom cases and with less errors compared to equivalent loads method [10]. 

The misalignment at the coupling has been studied along with the unbalance effect to 
show the synchronous forces and moments exists which depends on the type and amount of 
misalignment between the rotors [11]. 

The study of overhung rotors has gained importance, since the overhung rotor has more 
effect of its deflection due to its weight on its dynamic analysis. The gravity has a softening 
effect on the analysis of overhung rotors bringing the system towards the linearity [12]. The 
sudden unbalance effect like sudden mass loss in turbines can cause an impact effect which 
can excite the rotor to its natural frequency. Bode plots have been successful to analyze any 
sudden change in the unbalance of the rotors [13]. 


2. MATHEMATICAL MODELING 

The general equation of motion of overhung rotor system is given as: 
M q +C q QG q +Kq=0 


M q +(C+QG) q +Kq=0 

For calculation of Eigen values and Eigen vectors 
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where, M is the mass matrix, C is damping matrix, G is gyroscopic matrix, K is the 
stiffness matrix of the rotor system and q consists of the displacement and rotation terms at 
the nodes. The Gyroscopic forces and bearing properties varies with the shaft speed Q, hence 
matrix A and B will also varies as the shaft speed vary. For Timoshenko beam element in 
which shear and rotary inertia effects are included Fig.l, the element has two degree of 
freedom at each node which makes the element consisting of four degree of freedom. 
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Figure 1 Local coordinates for beam element 



Figure 2 Model representing unbalance and eccentricity respectively 

The unbalance response of the rotor system can be considered equal to the out of balance 
force response due to eccentricity. Hence for solving the unbalance response the forces 
generated due to eccentricity is considered for generation of equation of motion [14].The 
eccentricity of the rotor is shown in Fig.2 which can be compared with the unbalance of m u 
provided at the radius ofr 0 . The out of balance force due to eccentricity is mcQ . The lateral 
acceleration of m u is given as m u r 0 Q 2 . If out of balance force due to eccentricity and mass 
unbalance are equal, then 

m u r 0 = ms (5) 

where m is the total mass of the rotor and s is the eccentricity, for equality to exist 
between the above models s is small in comparison to r 0 and m u is small in comparison to m. 
Practically it is difficult to identify that whether the out of balance force exists due to 
eccentricity or mass unbalance. Study of unbalance is considered on Jeffcott rotor model with 
offset centre of mass from the shaft center line and the boundary conditions are applied for the 
overhung rotors. 

The equations of Jeffcott rotor in Fig.3 for free vibration supported by isotropic bearing 
are given as [14] 
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m x+ C xt x+ C xc y g + k xt x + k xc y 9 = 0 
m y+ C y , y~ C yc xe + k yt y - k yc xe = 0 
l d x e +l p n y e - Cyc y + C yr Xe — k yc y+ k yr x 0 = 0 
h y e - l P ^ x e + C xc x+ c xr y e + k xc x+ k xr y 9 = 0 

y 



( 6 ) 



Figure 3 Jeffcott rotor view perpendicular to x-axis and y axis respectively 

Where, x e and y e are displacement about x and y axis respectively. Due to the 

eccentricity s the force is induced in these free vibration equations due to which there are 
changes in the dynamics behavior of the rotor system which are given in equation (7). The 
variable notation is provided in table 1 . 

mu+ Ctu+ C c y e +k t u + k c y e = msCl 2 cos kit 

m v+C t y-C c xe + k t y-k c xo = msfl 2 sin Qt 

\ d x e + l p Cly e -C c v + C r xe-k c v + k r x e =0 

I d ye~^ p C^ x e + C c u+C r y e +k c u+k r y e = 0 


( 7 ) 


Table 1 Variables used in equation (7) 


Variable 

Specification 

Variable 

Specification 

k xl 

Stiffness of I bearing in x direction 


Damping of II bearing in x direction 

k x2 

Stiffness of II bearing in x direction 

Si 

Damping of I bearing in y direction 

k y i 

Stiffness of I bearing in y direction 

S2 

Damping of II bearing in y direction 

k y2 

Stiffness of II bearing in y direction 

C xt 

S 1 +S 2 

K 

k xl +k x2 

C * 

C yl+ C y2 

Kc 

-ak xl +bk x2 

C xc 

-ac xl +bc x2 

K 

a 2 k xl +b 2 k xl 

Cyc 

-ac yl +bc y2 

Kc 

-ak yl +bk y2 

C xr 

a C C x2 

K 

a 2 k yl +b 2 k yl 

Cyc 

a 2 c yl +b 2 c yl 

K 

k y ,+k y2 

K 

Diametral moment of inertia of the disc 

C xl 

Damping of I bearing in x direction 

K 

Polar moment of inertia of the rotor about z- 
axis 
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2.1. Unbalance Jeffcott Rotor Model Supported On Isotropic Bearings 

For research purpose the rotor is assumed to be supported on isotropic bearing i.e. the 
stiffness of the bearing in the each direction is same, the equations (7) can be modified as 

m u+ C t u+ C c y e + k t u +k c y e = mskl 2 cos kit 

m v+ C t v- C c xg+ k t v -k c xe = mskl 2 sin kit 

^d x e + ^pQye~ C c v+C r xe~k c v + k r x e = 0 

h ye- l P Q x e+ c c u+c r ye + K u+ K y d =° 


where, ^ xt= k v =^ t , 


( 8 ) 

Kr=K=K, K=K=K, C xt =c yt = c t , c„=c„=c, 


yr 


C xc = C yc = C c for isotropic bearing. In matrix form of equation (1): 
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2.2. Unbalance Jeffcott Rotor Model Supported On Anisotropic Bearings 

Practically the bearings on which the rotor is supported have different stiffness in every plane. 
So it is necessary to analyze the response for anisotropic bearings. The matrices for equation 
(1) are given as: 
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3. RESPONSE ANALYSIS AND RESULTS 

The overhung rotor system has coupling between the displacements and rotations. Letting 
r = u + jv and (p = y e ~ jx d equation (8) can be combined to form a single pair of equations 
for isotropic support. 
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m r+ c t r+ c c cp+ k t r + k c (p - msCl 2 e 


jClt 


J d <P~ <P + C c r+ C r (P + K r + K<P = 0 

The response of the system to out of balance force is given as [14] 

meD. 2 
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( 11 ) 
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Where D ~ d _ ^ + J c r ^ + X - m ^ 2 + J c t ^ ) _ 0 C C ^ + k c f 


Fig.4 shows the 3 dimensional view of the rotor system with two discs. The finite element 
modeling of the rotor-bearing-disc system has been done. This system is divided into 14 
elements as shown in Fig. 5. Element size between bearings is kept 54.6 mm and after node 
11 element size is 44.45 mm. The system carries two discs at node 6 and node 15. Shaft is 
supported by bearing at node 1 and node 11 representing an overhung rotor system. The 
unbalance mass m u is attached to disc at node 15 at a distance of r 0 from disc centre. Physical 
dimensions used for analysis are given in table 2. These parameters corresponds to the 
Machinery Fault system - Rotor Dynamics system (MFS-RDS) supplied by Spectra Quest. 
The material used for shaft is TGP steel and Aluminum for discs. For isotropic supports the 
value of k x is used as bearing stiffness in both directions and for anisotropic supports k x and 
k y are used as bearing stiffness. The analysis is performed through MATLAB programming. 



Figure 4 3-D model of the rotor system 
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Figure 5 FEM modeling of rotor system 
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Table 2 Physical parameters of the rotor system 


Parameter 

Value 

Parameter 

Value 

Total length of shaft 

723.9 mm 

k x 

7 x 10 7 N/m 

Distance between bearings 

546.1 mm 

kj 

5 x 10 7 N/m 

Diameter of discs 

152.4 mm 

O 

X 

II 

O 

■c 

5x 10 2 Ns/m 

Diameter of shaft 

19.05 mm 

m u 

5.43 gm 

Disc thickness 

14.27 mm 

r 0 

69.85 mm 


3.1. Isotropic Support 

Isotropic support has same stiffness in every direction. The response of unbalance is shown in 
Fig.7 through Bode plot diagram. In Fig.6 (a) we can see the response magnitude and phase 
change w.r.t. rotor spin speed for node 6 and node 15. The response of the unbalance which is 
provided at node 15 can be observed from the sudden fall in the response magnitude at 6860 
rpm where the phase also gets changed by 180 degree which corresponds to the same 
response in [13]. The other two peaks of 3960 rpm and 8320 rpm show response for critical 
speed as seen from the Campbell diagram (Fig.7 a) and they are for forward whirl. The 
response in x and y direction at the same node for isotropic supports are always out of phase 
as seen from Fig.6 (b). The first four mode shapes for the rotor system supported on isotropic 
supports is shown in Fig. 7 (b), it can be seen that odd number modes are backward and even 
number modes are forward. 




(a) Unbalance response for node 6 and 1 5 




(b) unbalance response for node 15 


Figure 6 Unbalance response for isotropic support 
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Forward whirl is red; Backward whirl is green. 



Nat Freq = 57. 8607Hz 



Nat Freq = 132.7309Hz 



Nat Freq = 61 4532Hz 



Nat Freq = 136.3346Hz 



(a) Campbell Diagram (b) Mode shape 

Figure 7 Campbell diagram and mode shape for isotropic support 


3.2. Anisotropic Support 

The response for anisotropic support is shown in Fig. 8. The two peaks in response magnitude 
at 3690 rpm and 8310 Fig.8 (a) are of critical speed as seen from Campbell diagram in Fig.9 
(a). There is sudden fall in response magnitude and change in phase at 6840 rpm which is due 
to unbalance response. There is some disturbance in response magnitude around 7830 rpm 
which can due to the mixed modes. The responses at 6840 rpm and 7830 rpm are of backward 
whirl which are shaded by yellow color in Fig.8 (b). The response for the same node in x and 
y direction is always out of phase as seen in Fig.8 (b). The odd number modes are backward 
and even number modes are forward as seen from Fig.9 (b). 




Yellow shading denotes backward whirl 



Rotor spin speed (rev/min) 


£ -too 


1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 
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(a) Unbalance response for node 6 and 15 


(b) Unbalance response for node 15 
in x and y direction 


Figure 8 Unbalance response for anisotropic support 
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Forward whirl is red; Backward whirl is green; Mixed whirl is black. 



(a) Campbell diagram 


Nat Freq = 54.6957Hz Nat Freq = 64.3632Hz 




Figure 9 Campbell Diagram and mode shape for anisotropic support 


4. CONCLUSIONS 

The response to unbalance of overhung rotor is studied on isotropic and anisotropic supports. 
The Bode plot is found to be useful in studying the unbalance response. There is sudden drop 
in the response magnitude due to unbalance for both the supports. Phase change is also seen 
where the amplitude suddenly drops due to unbalance in the rotor system. The response in 
anisotropic supports has mixed whirl for some range of rpm. The odd mode shapes are at 
slightly higher frequency for isotropic supports i.e. for backward whirl, where as the even 
modes shapes are at higher frequency for anisotropic supports i.e. for forward whirl. 
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